Pure Phase Decoherence in a Ring Geometry 
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We study the dynamics of pure phase decoherence for a particle hopping around an A'^-site ring, 
coupled both to a spin bath and to an Aharonov-Bohm flux which threads the ring. Analytic results 
are found for the dynamics of the influence functional and of the reduced density matrix of the 
particle, both for initial single wave-packet states, and for states split initially into 2 separate wave- 
packets moving at different velocities. We also give results for the dynamics of the current as a 
function of time. 
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I. INTRODUCTION 

The dynamics of phase decoherence is central to our 
understanding of those physical systems whose properties 
depend on interference. This is particularly evident when 
particles are forced to propagate around closed paths; 
phase coherence then makes all physical properties de- 
pend on the topology of these pathsii For this reason 
the quantum dynamics of particles on rings has been 
extremely important in our understanding of quantum 
phase coherence. Examples at the microscopic level in- 
clude the energetics and response to magnetic fields of 
molecules,^ as well as charge transfer dynamics in a vast 
array of solid-state and biochemical systems. There is 
evidence now for coherent transport around ring struc- 
tures even in some large biomolecules."^ At the nanoscopic 
and mesoscopic scale many ring-like structures, both con- 
ducting and superconducting,— show coherent transport 
around the rings, along with interesting Aharonov-Bohm 
style interference phenomena. We also note the impor- 
tance of closed loop structures in quantum information 
processing.— 

The interference around loops in all of these systems is 
very sensitive to phase decoherence. Questions about the 
mechanisms and dynamics of this decoherence are sub- 
tle, and have led to major controversies, notably in the 
discussion of mesoscopic conductors.^ A quantitative un- 
derstanding of decoherence processes in metallic systems 
and in superconducting "qubits" has yet to be attained 
(in both cases local defect modes clearly make the ma- 
jor contribution to phase decoherence at low temperature 
T) These controversies are examples of a wider prob- 
lem: typically in solid-state systems, low T decoherence 
rates are far higher in experiments than theoretical esti- 
mates based on the dissipation rates in these systems. 

These problems are complex because both decoher- 
ence and dissipation rates depend strongly on which 
environmental modes are causing the decoherence 
Delocalized modes (electrons, phonons, photons, spin 



waves, etc.) can typically be modeled as "oscillator 
bath" modesJ^"— In such models, decoherence goes 
hand-in-hand with dissipatiouj^ii^ in accordance with 
the fluctuation-dissipation theorem. However localized 
modes (defects, dislocations, dangling bonds, nuclear and 
paramagnetic impurity spins, etc.), which can be mapped 
to a "spin bath" representation of the environment ^li^ 
behave quite differently; indeed they often give decoher- 
ence with almost no dissipation. This is because although 
their low characteristic energy scale means they can cause 
little dissipation, nevertheless their phase dynamics can 
be strongly affected when they couple to some collec- 
tive coordinate - this then causes strong decoherence in 
the dynamics of this coordinate>i2ii£ The fluctuation- 
dissipation theorem is then not obeyed,^ and often these 
localized modes are rather far from equilibrium. 





FIG. 1: At left. An 8-site ring with nearest-neighbour hopping 
between sites. At right a potential U (R) with 8 potential wells 
(shown here symmetric under rotations by 7r/4), depicted as 
a contour map (with lower potential shown darker). When 
truncated to the 8 lowest eigenstates, this is equivalent to the 
8-site model. 

To understand how such non-dissipative decoherence 
processes work, it is then useful to look at models in 
which the environment causes pure phase decoherence, 
with no dissipation. As noted above, such models be- 
come particularly interesting when the decoherence is 
acting on systems propagating in 'closed loops'. Mod- 
els of rings coupled to oscillator baths have already been 
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studied.^ However such models, in which decoherence 
is inextricably linked to dissipation, do not capture the 
largely non-dissipative decoherence processes that domi- 
nate many solids at low T. On the the other hand pure 
phase decoherence has been studied in many papers, ^^^^"^ 
but not, as far as we know, the rather unique phenomena 
occurring on a ring. 

In this paper we study a model which embodies in a 
simple way both the 'closed path' propagation which is 
generic to quantum interference processes, and which in- 
volves pure phase decoherence coming from a spin bath. 
The model describes a particle propagating around a ring 
of A'' discrete sites, while coupled to a spin bath; we as- 
sume hopping between nearest neighbors. The model 
becomes particularly interesting if we also have a flux $ 
threading the ring (see Fig. [T]). The spin bath variables 
are assumed to be Two-Level Systems (TLS); these are 
ubiquitous in solid-state systems, and are the main cause 
of decoherence at low T in these systems. 

One can also study the problem of a continuous ring, 
but the discrete model is simpler, and is easily re- 
lated to diverse problems like quantum walks with phase 
decoherence or the dynamics of electrons in rings of 
quantum dotsi^ The Hamiltonian we will study has the 
general form 

= ^[A4cje^(^-+5^'=<-'^^)+i/.c.] (1) 

<ij> 

The operator c| creates a particle at site j; we assume a 
single particle only. The phase factors {A°j} result from 
the flux <I> threading the ring. In writing ([T]), we have 
assumed a symmetric ring, with A'^ sites, and assumed 
that the hopping matrix elements tij between sites i and 
j have simplified to a nearest-neighbour amplitude Aq 
(here X]<y> denotes a sum over nearest neighbours). 
This also means we can ignore any diagonal site ener- 
gies, since symmetry under rotations by angles 2-k/N 
means these energies are all the same. The spin bath 
variables {<Tfe} are Pauli spin-1/2 operators for the TLS, 
with k = l,2,....iVs. We emphasize immediately that 
these bath spins are, in real situations, often not spins, 
but instead the 2 lowest levels of localized modes in a 
solid (for example, as noted above, they could be defects 
or dangling bonds). 

The paper is organized as follows. In section II we dis- 
cuss the derivation of model Hamiltonians like ([T]) from 
more microscopic models, and the approximations which 
allow us to drop other terms that can also appear in the 
coupling of a ring particle to a spin bath. In section III we 
discuss the dynamics of the particle in the absence of the 
bath - this establishes a number of useful mathematical 
results. In section IV wc show how the dynamics of the 
reduced density matrix for the particle is derived in the 
presence of the bath, and give some results for this dy- 
namics. In section V we analyze the dynamics of a pair of 
interfering wave-packets moving around the ring, show- 
ing how pure phase decoherence destroys the interference 



between them. Finally, in section VI, we summarize our 
conclusions - since some of the calculations are quite ex- 
tensive, readers may want to look first at this section for 
a guide to the main results. The more technical details 
of the derivations in sections III and IV are given in an 
Appendix. 

II. DERIVATION OF MODEL 

Consider first an A^-site ring system without a bath. 
In site representation, this typically has a "bare ring" 
model Hamiltonian 

Ho= h'^Ic, e^^^ + H.c] + ^ e,c]c, (2) 

<ij> 3 

This " 1-band" Hamiltonian is the result of truncating, to 
low energies, a high-energy Hamiltonian of form: 

Hv = ^^{V-A{n)f + U{R) (3) 

where a particle of mass M moves in a potential U{R) 
characterized by A^ potential wells in a ring array (see 
again Fig. [1]). Then ej is the energy of the lowest 
state in the j-th well, and tij is the tunneling ampli- 
tude between the i-th and j-th wells (which we take here 
to be nearest neighbours). In path integral language, 
this tunneling is over a semiclassical "instanton" trajec- 
tory Rins (r) , occurring over a timescale ~ 1 /i^o (the 
"bounce time")j^ Here fio (the "bounce frequency") is 
roughly the small oscillation frequency of the particle in 
the potential wells. In a semiclassical calculation, the 
phase A"j is that incurred along the semiclassical tra- 
jectory by the particle, moving in the gauge field A(R). 
For a symmetric ring the site energy Sj — > Eq, Vj, and we 
henceforth ignore it. 

Consider now what happens when we couple the par- 
ticle to a spin bath. The spin bath itself, independent of 
the ring particle, has the Hamiltonian 

HsB = Y.^k-<T, + V^,P,aZ4 (4) 

in which each TLS has some local field acting on 
it, and the interactions V^^, are typically rather small 
because the TLS represent localized modes in the envi- 
ronment. The most general coupling between the ring 
particle and the bath has the form 

j 

+ Y,{G%{au)c\c,+H.c.)] (5) 

in which both the diagonal coupling F^^ and the non- 
diagonal coupling G*j are vectors in the Hilbert space of 
the k-tYi bath spin. We shall see below, when considering 
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the origin of these terms from microscopic models, that 
very often we can write the total Hamiltonian as 

(6) 



where H^and = Ho 



H, 



bajid 



H — Hhand + HsB 

Hint takes the form 



H.c] 



(7) 



in which the diagonal couplings to the spin bath assume a 
"Zeeman" form, of strength |7^|, linear in the {crk}, and 
the non-diagonal couplings appear in the form of extra 
phase factors in the hopping amplitude between sites. 

Before we consider the microscopic origins of this 
model, let us note how it simplifies when we assume the 
symmetry under rotations by 2tt/N noted above (so that 
the site energy ej is dropped, and tij — >■ A^, with nearest- 
neighbour hopping only). It is then natural to write A° 
as 



4° 



-H ■ X R^. = <i>/iV 



(8) 



for j = i + 1 (we now use MKS units, and put h = 1). 
Here, H is the magnetic field, and R^ is the radius- vector 
to the ith site; in cylindrical coordinates 



Rj 



iRo,e,) 

2nj/N 



(9) 



for a ring of radius Rq. Fourier transforming from the site 
basis to a momentum basis for the couplings, we define 
quasi-momenta A:„ = 2Tm/N, with n = 0, 1, 2, A'^ — 1, 
for the particle on the ring, and define operators 



-E 



N E 



N 

fc„ = ^ , " = 0,1,...,A^-1 . 
We can write the free particle Hamiltonian as 



(10) 



2Ao^cos(fc„-(f>/iV)ct^Cfc^^ 



(11) 



Then in this basis we can write: 

V.nt =Y.T. [Fn{^k)p{K) + Gi{au)cljuJ\ (12) 

k n 

where p(/c„) — X^n' +k density operator 

in momentum space for the particle, and the Fourier- 
transformed interaction functions are 



e 



(13) 



In this basis the band Hamiltonian Hhand has a dispersion 
which is a functional of the bath spin distribution: 



Hband = E E'^fc.. [""fcl'^L ^fc^ 



k n 



(14) 



and in which the 'band energy' ek„ [""fe] and the 'scatter- 
ing potential' w„[crj.] are now both functionals over the 
spin bath coordinates {cfc}: 



Efe, 



VnlCTk 



E^«( 



(15) 



Under many circumstances one can assume that this 
symmetry under rotations also applies to the bath cou- 
plings, so these no longer depend on site variables, ie.. 



F^, and G% G'' . The results then simplify a 



great deal; Gn{(Tk) 2G''{(Tk) cosfc„, and F^ . 

Now let us consider the microscopic origin of this model 
(ie., before truncation to the lowest band). The most 
obvious interaction between the particle moving around 
the ring and a set of bath spins has the local formi^: 

fc 

^ 5]ifL(R,'Tfc) (16) 

fc 

where F(r) is some vector function, and is the position 
at the fc-th bath spin. The diagonal coupling Fj, or its 

linearized form 'y^, is then easily obtained from (jl6p when 
we truncate to the single band form. But the term (jl6p 
must also generate a non-diagonal term, which is more 
subtle. We can see this by defining the operator 



exp [—i/h 



(17) 



where the particle is assumed to start in the i-th potential 
well centered at position R^, at the initial time t™, and 
finish at position Rj in the adjacent j-th well at time 
r/; the intervening trajectory is the instanton trajectory 
(which in general is modified somewhat by the coupling 
to the spin bath). Now we operate on CTfc with Tf'j, to get 



rj-^k 



•k I 



'k I 



(18) 



where we note that both the phase (^^^ multiplying the 
unit Pauli matrix cr° , and the vector ol^ multiplying the 
other 3 PauH matrices al,a^,a^, are in general com- 
plex. In this way the instanton trajectory of the par- 
ticle acts as an operator in the Hilbert space of the A:-th 
bath spiuiiSiS! Note that one important implication of 
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this derivation is that typically \a^\ ^ 1, in fact expo- 
nentially small, since the interaction energy scale set by 
\F{R— rk)\ is usually much smaller than the "bounce 
energy" scale Mlo set by the potential C/(R), ie., the 
tunneling of the particle between wells is a sudden per- 
turbation on the bath spinsj^ Detailed calculations in 
specific casesiSi^^i^ show that ja^:' | ~ 7r|a;[:' |/2rio in this 
'sudden' regime, where u/j! — — 7)^, is the change in 
the diagonal coupling acting between the particle and the 
fc-th bath spin when the particle hops from site i to site 
i (this result can be found directly from time-dependent 
perturbation theory in the sudden approximation). 

From these considerations we see that, starting from a 
ring with the particle-bath interaction given in (I16p . we 
will end up with an effective Hamiltonian for the lowest 
band of the form given in ([7]) , in which the non-diagonal 
interaction G^{(rk) in ([5]) has assumed a rather special 
form. 

One can in fact have a more general form for G^j{crk) 
in the lowest-band approximation, provided one also in- 
troduces in the microscopic Hamiltonian a coupling 



(19) 



to the momentum of the particle. This can include var- 
ious terms, including functions of P x cr^ and P ■ CTk] & 
detailed analysis is fairly lengthy. The main new effect of 
these is to generate terms in the band Hamiltonian which 
couple the spins to the amplitude of tij as well as to its 
phase; these do not appear in ([7]). 

In any case, if we know C/(R), F{R — r^), and 
G{P,(Tk), we can clearly then calculate all the param- 
eters in the generic model Hamiltonian, using various 
methodsj ^°'^^ However we are not interested here in the 
generic case, since our main object is to study the dy- 
namics of decoherence in a ring model which contains 
only phase decoherence. We therefore make the follow- 
ing approximations: 

(i) We drop the interaction V^/^ , between bath spins 
(often a very good approximation, since interactions be- 
tween defects or nuclear spins are often very weak), and 
also neglect the local fields acting on the {(t^}. Thus 
we make Hsb = 0. 

(ii) We drop the momentum coupling G{P, cTk) en- 
tirely, and in the band Hamiltonian ([7]) we drop the di- 
agonal interaction 7^. This implies that the energy of 
the A:-th bath spin does not depend on whether the j-th 
site is occupied. We make this approximation (in many 
cases not physically reasonable) only because we wish 
to study phase decoherence without the complication of 
energy relaxation. 

(iii) We assume a symmetric ring, so that Sj — )• and 
tij — >• Ao as before; and we absorb the phases (j)]^ — >■ 4'k 
into a renormalization of Aq (from X^fc 'f'k): ^-nd of 
A° (fromEfeRe^fe). 

The resulting model i/^ is then just that given in ([Ij. 
This turns out to be explicitly solvable, and reveals some 



important properties of phase decoherence. We will usu- 
ally assume the parameters a'j^ are small, in line with the 
remarks above (although the net effect of all of them may 
be very large), and we will also usually specialize to the 
case a]^ — >■ cc^, consistent with a completely symmetric 
ring. 

Finally, let us briefly compare with the kind of Hamil- 
tonian one would expect for a particle on a ring coupled 
to an oscillator bath. Let us assume a set of oscillators 
with Hamiltonian Ho + Hose + Hint, where Ho is again 
the free particle hopping Hamiltonian, coupling to a set 
of No oscillators with Hamiltonian 



^° 1 p2 



(20) 



9=1 



In general there will be diagonal couplings {Vj{q)} 
and non-diagonal couplings {Uij{q)} between particle 
and oscillators. We could also have a coupling to 
the oscillator momenta - however in this case one can 
make a canonical transformation^® which transforms 
this back into a coupling to the {xq}. Typically the 

— 1/2 

couplings {Vj{q),Uij{q)} ~ 0{No ). We note here 
that in many microscopic models of this kind, the 
coupHngs {Vj{q),Uij{q)} are actually also strong func- 
tions of temperature, either because the underlying ef- 
fective Hamiltonian is strongly T-dependent (eg., in a 
superconductor) j-^ or because the coupling to the oscil- 
lators is non- linear (eg., in the coupling to a soliton)4S 
If we restrict the problem to rotationally invariant cou- 
plings on the ring, then we can write 



Hint — 



E 



E 

<ij> 



iu4c,+H.c.) + J2Vqc]c 



Xq (21) 



where Uij{q) Uq, Vj{q) -> Vq, and the sum J2<ij> is 
over nearest neighbours. It is then straightforward to go 
through the same manipulations as in (IT2|) - (IT5t . to get a 
renormalised band which is a functional of the {xq}. 

In these results there is no connection between the ring 
sites and the space in which the oscillators are supposed 
to exist. However in many cases the oscillator displace- 
ment field Xj can be defined at each site j of the ring; 
the coupling then reduces to 

q j 

= E"(f^J-^J')S^^J-^(^J') (22) 
jj' 

in which Rj, Tji are site vectors on the ring, and Xq is 
now the Fourier transform of Xi. 



III. 



FREE BAND PARTICLE DYNAMICS 



We first consider the dynamics of a free particle in 
some initial state moving on the symmetric A^-site ring 
described by Ho in (fTT|) . with no bath. 
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For this free particle the dynamics is entirely described 
in terms of the bare 1-particle Green fmiction 

^ ]_'^(,-i2Agtcos{k^-'i>/N)^ikJj'-j) ^ f^23) 
n 

which gives the amplitude for the particle to propagate 
from site j' at time zero to site j at time t. These paths 
are rather simple (see Fig. [5]); they can be labelled by 
the initial and final sites, and by the winding number of 
the path around the ring. 




FIG. 2: Colour online: A particular path in a path integral 
for the particle, shown here for an A'^ = 3 ring. This path, 
from site to site 1, has winding number p — 1. 

The 1-particle Green function can be evaluated in var- 
ious ways (see Appendix) ; the result can be usefully writ- 
ten as 

P— — 00 

(24) 

where is a sum over winding numbers. The "return 
amplitude" Goo(i) is then given by 

p 

= Y,e~''"^lNp{-2i\t) (25) 
p 

where in the last form we use the hyperbolic Bessel func- 
tion. 



where p^p' are the winding numbers (see Appendix for 
the derivation for a general initial density matrix). This 
form has a simple physical interpretation - the particle 
propagates along pairs of paths in the density matrix, one 
finishing at site j and the other at site j', and the order 
of each Bessel function simply gives the total number of 



It is often more useful to have expressions for the den- 
sity matrix; even though these depend trivially for a free 
particle on the Green function, they are essential when 
we come to compare with the reduced density matrix for 
the particle coupled to the bath. One has, for the 'bare' 
density matrix operator of the system at time 

p°{t) = e-^^°*p°(0)e'^°*. (26) 

Thus, suppose we have an initial density matrix pi™'' = 
{l\p{t = Q)\l') (where / and V are site indices), then at a 
later time t we have 

= pPG^t)G''^,At)\ (27) 
where we use the Einstein summation convention (sum- 
ming over 1,1'). This equation defines the propagator 
K°j, ii{t) for the free particle density matrix, as 

i^°,.,,„,(i) = G°,(t)G°,,,(i)t. (28) 

In the main text of this paper we will almost always 
quote results for the special case where the particle be- 
gins at t = on site 0. In the case of the free particle, 
this means that = 5qi5i'q, and only the propagator 
matrix K°y ^^{t) enters the results; then we have 

{j\p°m') ^ if°-,,oo(i) = G°,{t)G%{t)\ (29) 

In the Appendix we give the results for an arbitrary ini- 
tial density matrix. 

The evaluation of the time-dependent density matrix 
for the free particle turns out to be quite interesting 
mathematically. As discussed in the Appendix, one can 
evaluate p"^, (t) as a sum over pairs of paths in a path 
integral, to give a double sum over winding numbers, or 
else as a single sum over winding numbers. Consider first 
the double sum form; again, for the special case where 

pI™'' = <5o;<5;'o i^^'^ particle starts at the origin), this can 
be written as 



(30) 

I 

sites traversed in each path, with appropriate Aharonov- 
Bohm phase multipliers for each path. 

If one instead writes the answer as a single sum over 
winding numbers, again assuming = <^o;'^;'0' S^^' 



J 



pp' 
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m— p' — — oo 



where as before the {km} are the momenta of the par- 
ticle eigenfunctions. The physical interpretation of this 
form is less obvious, but the sums are much easier to 
evaluate since they only contain single Bessel functions 
instead of pairs of them. Thus wherever possible we re- 
duce double sum forms to single sums. Notice that for 
these finite rings, the bare density matrix is of course 
strictly periodic in time. Notice also that the diagonal 
elements of p{t) are generally periodic in $. However, the 
off-diagonal elements are only periodic in ^/N. In con- 
trast, e**^-'"-' {j\pit)\j') is periodic in $, with period 
2tt. This latter is the quantity needed for calculating the 
currents, as we will see below. 

From either Gjj'{t) or Pjjr{t) we may immediately 
compute two useful physical quantities. First, the prob- 
ability Pjoit) to find the particle at time t at site j, as- 
suming it starts at the origin; and second, the current 
Ijj+i{t) between adjacent sites as a function of time. 

Looking first at the probability /'°o(0' 

pm - m^m) - i^oWP (32) 

which from above can be written in double sum form as 

Pjoi^) =E JNp+j{'2^ot)JNp'+j('2^ot) 
PP' 

X e-*^(p'-f)(*/^+f) (33) 
or in single sum form as 

N~l 00 

P^oit) E E e-(*+^'='"/^) 

m— P— — 00 

X Jjvp[4A„tsin(A:„/2)] . (34) 

One may also compute moments of these probabilities. 
These are not terribly meaningful for a small ring, be- 
cause any wave-packet will be spread around the ring. 
However for a large ring they can be useful- for example, 
the 2nd moment J^j i^-fjo(0 tells us the rate at which an 
initial density matrix spreads in time, provided the spa- 
tial extent of the density matrix is much smaller than 
the ring circumference. Coherent dynamics will then 
manifest itself as ballistic propagation of an initial wave- 
packet. 



From these general expressions it is hard to see what 
is going on. To give some idea of how the probability 
density behaves, it is useful to then look at these results 
for a small 3-site ring, where the oscillation periods are 
quite short. One then has, for the case where the particle 
starts at the origin, that 

Pjoit) = ^(1 + (33^,0 - l)[Jo(2AoV3t) 

OO 

-t- 2 ^ J6p(2Ao Vst) cos(2p$)] 

p=i 

OO 

+ (<5^- 1 - <5,, 2)2^3 ^ J,^_,{2AoV3t) sin((2p - 1)$)) . 

(35) 

In Fig. |3]the return probability P§Q{t) is plotted for the 
case = 3, using ([55)1 . From the results one striking 
feature immediately emerges - we see that the periodic 
behaviour depends strongly on the flux $. This flux de- 
pendence illustrates the way in which the flux controls 
the particle dynamics, by acting directly on the particle 
phase. In section V we will see how this also happens 
when one looks at interference between 2 wave-packets; 
and in sections IV and V we will see how decoherence 
washes out the flux dependence of the particle dynamics. 
Thus the flux dependence of the particle dynamics very 
effectively measures how coherent its dynamics may be. 

Turning now to the current 1° j+i {t) from site j and site 
J -|- 1, this is given from elementary quantum mechanics 
by 

Il^^S) = 2Im[A„e-*/^p°,.+i(t)] 

= zA„(e'*/^p°+i,^.(t)-e-*/^p°^.+i(i)) (36) 

where the flux per link appears in each contribution. 
Again, one can write this expression as either a double 
sum over pairs of winding numbers, or as a single sum 
(see Appendix for the general results and derivation). For 
the case where the particle starts from the origin, these 
expressions reduce to 
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(2Aot)cos[(-7V + <I>)(p'-p)] 

OA h 

m— p 

I 



for the double and single sum forms respectively. 



Pj 




FIG. 3: Colour online: Results for the free particle for = 3 
and for a particle initially on site 1. Left: The probabilities to 
occupy site 1 (full line) , 2 (large dashes) , and 3 (small dashes) . 
Right: the current from site 1 to site 2. Top: $ = 0. Bottom: 
$ = tt/2. 

Again, the currents across any links must be strictly 
periodic in time for this free particle system; and again, 
it is useful to show the results for a 3-site system. For 
this case = 3, and assuming that the particle begins 
at the origin, we find 

2 

IS,i - ^Im^^J3p+i(4Aotsin— ) 

m=i P (38) 

X e— (3p+i)/3i3p+i cos[(y + $)p] 

which we can also write in the form 

IS,i = ^Im^ J3p+i(2V3A<,t) cos[(^ + m 
p 

2 

^ -Sp+l ^ (g-i^(3p+l)/3 _^ g-427r(3p+l)/3^ ^gg^j 
m— 1 

Now let us write (e-*'^(3p+i)/3 + e"'2,r(3p+i)/3) ^ 
(— l)Pe~*'^/3 + e~^"/3. If p is even, this becomes —iV^ 
and cos[(^ + $)p] ^ (-l)3p/2pos($p); If p is odd, it 
becomes -1 and cos[(^ + <I>)p] = {-1)^'-p-^'^/'^ sm{^p). 



Therefore, we have 
2 °° 

/o°i=-A„ ^ J3p+i(2A„V3t)if(p,<I') , 

p— — oo 

K{p, $) sin(|5$) if p = odd , 

Ar(p, $) = \/3cos(p"I>) if p = even . (40) 

These results are also shown in Fig. [3] Notice that in 
this special case the result is periodic in $; this is not 
however true for a general initial density matrix p^r\ 
when the result is periodic in $/iV. 

IV. RING PLUS BATH: PHASE AVERAGING 

We now wish to solve for the dynamics of the particle 
once it is coupled to the bath, via the Hamiltonian ([T|). 
This is done in general by integrating out the bath spins, 
to produce expressions for the reduced density matrix 
of the particle. In this section we first show how this is 
done, and then give results for physical quantities (in par- 
ticular, the probability Pjo{t) and the current Jjj+i(i)). 
Finally, we briefly compare the results to the behaviour 
one expects for a ring coupled to an oscillator bath. 

A. General results 

As shown in the Appendix, the reduced density matrix 
for the particle obeys the equation of motion 

P.yW = E%'^«'W^i? (41) 

where JCjji ji' (t) is the propagator for the reduced density 
matrix. This latter can be written in the form of a double 
sum over winding numbers 

^Y.^"rM'ip,p';t)F^rip,p') (42) 
pp' 

where the function Kjj, ii,{p,p';t) is the free particle 
propagator for fixed winding numbers p,p' (so that 
K°,',wit) = Epp'K;;^',iAP,P'-^ty, see the Appendix, 
eqtn. (jASp et seq ). All effects from the spin bath are 
then contained in Fjj,{p,p'), which we will call the "in- 
fluence function" . The remarkable thing is that this func- 
tion depends only on the initial and flnal states, and on 
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the winding numbers - all other aspects of the two paths 
involved in the density matrix propagation have disap- 
peared. As explained in the appendix, this is a particular 
feature of the pure phase decoherence being treated here. 

The form the influence function takes depends on what 
kind of averaging we do over the bath. To discuss this, 
let us first discriminate between two different ways of 
averaging over the bath, as follows: 

(i) The first and most obvious case is where the a™" 
are considered to be a set of fixed couplings, for a specific 
single ring. In this case the average is only over the bath 
states; we will denote this bath average by < .... >. Often 
it will only involve a thermal average over the bath states. 

(ii) However it is often the case that one is either in- 
terested in an ensemble of rings, all having the same free 
particle Hamiltonian but with the cxk possibly varying 
from one ring to another, or a single ring in which the 
values of the couplings cxk are indeterminate. In this case 
it makes sense to define a probability distribution P{a) 
over a coupling variable a. One then must average not 
only over the bath states themselves, but also over the 
bath couplings. We will denote this double average by 

<< >>, to signify the average over both the bath 

states and the probability distribution; and the influence 
function for this case will be written as F^''j,{p,p'), with 
the bar over the F signifying that an average over cou- 
plings is being done as well. 

In general the results for the dynamics of the density 
matrix and the current, and their dependence on the in- 
fluence function, may be quite complicated. Thus, before 
we begin quoting results, it is useful to note what are the 
important parameters in the problem. We will only con- 
sider here the simplest completely symmetric case where 
a^" — > otk for all links {mn}; and we will assume that 
|Q!fe| <^ 1 for all k, as discussed in section II. Now in the 
previous literature for this case of pure phase decoher- 
ence, it has been usual to deflne a 'topological decoher- 
ence' parameter-i2i^ 

A = ^El"fcl' (43) 

k 

which provides a measure of the strength of the pure 
phase decoherence lis If the number Ng of bath spins is 
large, then we can have A ^ 1; this is the limit of strong 
phase decoherence. 

However we shall see in what follows that on a ring 
it is often more useful to define a parameter Fq{p) that 
also depends on a winding number p. The form of this 
parameter depends on which of the two bath averages is 
performed. In the case where only an average over the 
bath states is performed, we have 

FoiP) = l[cos{Np\ak\) (44) 
fe 

which defines a rather complicated function of the fixed 
bath couplings. The strong decoherence limit for this 
case is defined by the parameter A defined above. 



In the case where we also perform an average over the 
bath couplings, we have 

MP) = n / 'i^kPiak) cos{Np\ak\) (45) 

k 

The result then depends on what form one has for the 
distribution function P{a.k). In what follows we will use, 
as an example, a Gaussian distribution, given by 

P(|afe|) = e-l"^l'/2A„/y2^ (46) 

so that 

Fo(p) = e-^^'p"/^ X^NsX, (47) 

The limit A — >■ oo is the "strong decoherence" limit for 
this distribution, where we have Fo{p) — )■ SpQ. How- 
ever we will see below that it is convenient to think of 
the strong decoherence regime for the present problem 
as that for which the particle dynamics is independent 
of fiux - we will see that this happens already for quite 
small values of A. 

We can see why these functions enter by considering 
the forms for Fjj,{p,p') and Fj''j,{p,p') that enter into 
physical quantities. In the appendix the full expressions 
for these are derived; but here we will again only use 
them for the case where /cj™'' = SqiSi'o, ie., the particle 
starts at the origin, and so only the function Fjji(p,p') = 
-Fj^j] {p,p') comes in. We will also again assume the purely 

symmetric case where — )■ otk for every link. 

Let us first consider the case of fixed bath couplings. 
In this case the form of the influence function reduces to 
(see Appendix): 

Fjj,{p,p') = (e-»^[(P-p')+(J-j')12:.ai--'^i-) (48) 

Notice that Fjjr{p,p') is a function only of the distance 
j — j' between initial and final sites, and of the difference 
p ^ p — p' in winding numbers. Writing this now as 
Fjji (p) , let us evaluate it by assuming the usual thermal 
initial bath spin distribution. Since all the bath states 
are degenerate, then at any finite T all states are equally 
populated; we then get: 

Fjj'iP) = l[cos{{Np + j~f)\ak\) (49) 

k 

Other initial non-thermal distributions for the spin bath 
states are also easily evaluated from pSj) . 

B. Physical Quantities 

From expressions like (j49p one can now write down 
expectation values of physical quantities as a function of 
time. The simplest example is the probability for the 
particle to end up at some site after a time t, having 
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started at another. Thus, eg., the probabihty Pjo{t) to 
move to site j from the origin in time t is now given by 

= JNp+j{2Aot)JNp'^j{2Aot) 
pp' 

X e-'Nip'-pK<s>/N+i) Fo{p,p') (50) 

which is a simple generahzation of the free particle result 
in p3p : we note that only the term 

Fo{P,P') = llcos{N{p^p')\ak\) (51) 

k 

in the influence function survives in this expression. Since 
this function depends only on the difference p — p' , it 
is identical to the function ^o(p) defined in (H^ above 
(letting p — p). We shall see below that the ring current 
is also controlled by this same function. Note that it 
has a complex multiperiodicity, as a function of the Ns 
different parameters A'^pjccfcl; we do not have space here 
to examine the rich variety of behaviour found in the 
system dynamics as we vary these parameters. 

Now let us consider the case where we also average over 
the bath couplings. One then finds (see appendix) that 

Fn'iP) = n/ rfafcP(afc)(e-*^[^^+(-'-^')l"-'^'=) (52) 

k 

In the symmetric case we can treat each bath spin in the 
same way, and simply use a distribution function P(|a|), 
the same for all the different {(Tk}- Then we can treat 
everything in terms of this single average, over a single 
representative spin cr from the bath. Then, eg., for an 
initial thermal ensemble for the bath spins, this gives 

Foo'{p) = [ j dcxP{cx)cos{{Np + j-3')\oc\)f^ (53) 

To give something of the flavour of this case, we use 
the Gaussian distribution for the P(|q;|), given by (pS)) 
above. Then, for the thermal ensemble just given, we 
have 

(p) - exp[-A(iVp + j' - jf/2] (54) 

It is then immediately obvious that the result for the 
probability for the particle to go from site to site j in 
time t is the same expression as (I50p above, but now with 
Fq{p) instead of -Fb(p)- 

To see how this behaves, let us take the specific case 
where TV = 3 again. Then for this 3-site ring one has, for 
example, that 

AoW = J(l + 2[Jo(2Ao%/3t) 

OO 

+ 2^ J6p(2AoV3t)cos(2p$)Fo(6p)]) . 

(55) 




FIG. 4: Colour online: Plot of Pjo{t) for a 3-site ring, for 
a particle initially on site 1, in the intermediate decoherence 
limit, with A = .02. Left: The probability to occupy site 
(full hne), 1 (large dashes), and 2 (small dashes). Right: the 
current from site to site 1. Top: $ = 0. Bottom: $ — 71/2 . 



To analyse this result, note that for x ^ (6p)^, we can 
use Jqp{x) « {—iy^/2/{nx) cos{x — 7r/4). The function 
Fq{P) decays with p, and becomes negligible for large 
enough p. For example, Eqtn. (|47| implies that {6p) < 
e-^° for p > p^^^, with p^^^ = V5/9A. Neglecting 
these terms in the sum in Eq. (j55p we conclude that for 
2AoV3i > {QPmaxY have e.g. 

1 2yl 

Pioit) « o [1 + ; cos(2A„ V3t - 7r/4)] , 

OO 

yl = l + 2^(-l)Pcos(2p$)Po(6p) . (56) 

The sum in the amplitude A reduces to '^{—ly Fo{&p) , 
for $ = 0, and to Fo{Qp), for $ = 7r/2. Clearly, switch- 
ing from $ = to $ = 7r/2 causes a large increase in A. 
Notice that the inverse Fourier transform of the ampli- 
tude A{4>) can be used to measure the decoherence func- 
tion Fq{6p). Results for this low decoherence regime are 
shown in Fig. 21 

As A increases, p^ax decreases, and Eq. (1551) applies at 
shorter times. Remarkably, if A > 0.1 the whole sum be- 
comes negligible, and we have already reached the strong 
decoherence result where the result is $— independent. 
The result is shown in Fig. 13 Thus, if we define the 
'strong decoherence' regime as that where all results are 
flux-independent, then it is reached for very low values 
of A. We emphasize here that the detailed form of the 
results, as well as the decoherence strength required for 
flux-independent dynamics, depends strongly on the form 
we adopt for either Fq{p) or F{p); we do not have space 
to explore this question here. 

Turning now to the current through the ring, we gen- 
eralize the free particle results in the same way as above. 
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FIG. 5: Colour online: Plot of Pjo{t) for a 3-site ring, for 
a particle initially on site 1, in the strong decoherence limit. 
Left: The probability to occupy site (full line), 1 (large 
dashes), and 2 (small dashes). Right: the current from site 
to site 1 (compare Fig. [31). The results do not depend on $. 



Quite generally one has 



where we average the operator 



(58) 



{t) = i{Ajj+ipj+ij{t) - Aj+ijpjj+i{t)) (57) 



over bath states, with fixed bath couplings - the case 
where one also averages over an ensemble of bath cou- 
pUngs is a by now obvious generalization of this. This ex- 
pression is evaluated in detail in the Appendix; as noted 
there, the result is more complicated than it seems, be- 
cause the density matrix depends implicitly on both the 
initial state, and on the full details of the propagator for 
the density matrix. Here we consider only the special 
case where the particle starts from the origin, and the 
fully symmetric case a^j^ — > ct^. Then one has, for the 
case of a bath state average only, that 



it) 



2A 



N-l 



^ ^ ^ Jnp+1 (4Ao< sin ■ 

m— p 



(59) 



r 



with a similar result for the current Ij j^i{t) arising in 
the case where one also averages over bath couplings, 
with Fq{p) then replaced by Fo{p). One can also analyze 
this result as a function of time, and of the decoherence 
strength, the ring size, and the flux - there is no space 
for this here. To nevertheless give some flavour for the 
results, consider again the 3-site ring, for the coupling 
averaged case, in the strong decoherence limit. Current 
then only flows in regions where the initial density matrix 
is inhomogeneous; for some general initial density matrix 
one finds 



^A„(,<:;) 



(60) 



where p[l7^ is the initial density matrix. Again we see 
that the result is completely independent of the flux. 



C. Comparison with Oscillator Bath 

To gain some perspective on the results just given, it is 
useful to compare with what one might expect for a ring 
particle coupled to an oscillator bath. The differences 
are both formal and physical, and both are important. 
Here we simply sketch these - a more detailed study of 
this rather complex problem will appear elsewhere To 
specify the formal problem completely, one needs first to 
define 'spectral functions' for the couplings between the 
oscillator bath and the ring particle.t^ These couplings 
were defined earlier, in (j2ip : Fourier transforming them 



in the same way as we did for the spin bath couplings, 
we then define the spectral functions as: 



S{UJ ~ LOq) 



4(' 



W) =- 2^ — 6{uj~UJq) 



(61) 



In many cases the non-diagonal function Jpiuj) can be 
neglected compared to the diagonal Jp(w), and we will 
assume this here. Jp{'^) can take many forms; the 
most commonly analysed is the " Ohmic form" , where 
Jp(a;) = rju) at low frequency, but this form is very useful 
for systems coupled to an itinerant electron bath, it is 
inappropriate for insulating systems (where a more accu- 
rate low-w form is the "superOhmic" form Jp(w) ~ w*^, 
with fc > 1). In addition, there is often significant low- 
energy structure in Jp(w), not describable by a simple 

power-law form; and in many cases Jp(w) also depends 
strongly on temperature T. 

Defining the influence functional J-[Q, O'] in the usual 
way for general paths 0(i),0'(<) (cf. eqtn. (|A19p of the 
Appendix), we can write 



J"[e, 6'] = exp 



dti I dt2 X 



[ip{ti)Dp{ti - t2)^{t2) + iTp{ti - t2Mh)i;{ti)] (62) 
where we have defined the sum and difference angular 
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variables 



V'(t) = (e(t) + e'(t))/2 



The behaviour in time of Dp (t) can be quite complex, and 
varies strongly with the form of Jp(w), and with temper- 
ature; the details of this behaviour have been reviewed 
(63) extensively J^i^ 



and the oscillator propagator Vp^t) = Dp{t) + iTp{t), 
with 



^V^(l - coswi coth^ 



In the same way as for the spin bath, we may now 
construct expressions for the reduced density matrix, and 
physical correlation functions derived therefrom, by sum- 
ming over all paths; this is done in a simple general- 
isation of methods developed for the spin-bosoit^ and 
(64) Schmid^'^''^^ models. For example, the probability Pno{t) 
takes the form 



OO OO 



P„o(t) ^ (-l)'-"-^''e^*(P+"/^)A2' / dt 

p=-ool=\n+Np\ 



r 
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JO 



T.T.^iMd'^rhiU}) (65) 



written as a sum over winding numbers p and the num- 
ber of intersite hops I. In this expression the influence 
functional has now become a function F{{^r}, {Xr}] {U}) 
of the times ti at which the particle hops, and of two sets 
of 'charges' {g^} = ±1, Wr} = ±1. These charges are 
defined in terms of the sum and difference paths by 



/ 



(66) 



so that the describe hops in the 'centre of mass' part 
of the density matrix, and the <Tr are hops in the 'differ- 
ence' or off-diagonal elements of the density matrix. The 
general form of -F({Cr}, {Xr}; {U}) is 



21 



Fi{ir},{Xr};{tl}) = 6(2,2 -Y,'lr)SiY,<^r) 

r—1 r—1 

X exp — \D^(tr — tr')qrO'r' + iTitj. — tr')(7rO'r'] 



t' <r 



(67) 



and we get the well-known oscillator-mediated interac- 
tions between the charges, familiar from the spin-boson 
and Kondo problems. Thus from the formal point of 
view, for a ring particle coupled to either an oscillator or 
spin bath, the principal difference between the two cases 
is the existence, in the oscillator bath case, of retarded in- 
teractions between particle hops at different times, whose 
form depends on Jp(w) and on T. Just as in the spin- 
boson and Schmid models, the interactions between the 
charges in the Ohmic case eventually cause a zero tem- 
perature Kosterlitz-Thouless binding transition between 



the charges, which localizes the particle at one site in the 
ring. This happens at a critical Ohmic coupling strength 
r] ~> rjc — h/2iT, independent of the ring size*2i When 
rj < rjc, the particle dynamics is strongly diffusive - we 
do not go here into the details of how the dynamics varies 
with ?7, with N , and with temperature T . In the super- 
Ohmic case there is no localization transition, no matter 
how strong the coupling; the analysis of this case is very 
lengthy 

None of these features has any formal counterpart in 
the coupling to a spin bath. In the present case the spin 
bath results are entirely independent of T, because all 
bath levels are degenerate. Even when this is not the case 
(ie., when we add back the local fields {hfc}, so that the 
decoherence becomes temperature dependent), the only 
way that interactions can be generated between different 
bath spins is through their coupling to the particle itself 
- there is no analogue to the propagator 'Dp{t). 

The key physical difference between this ring-oscillator 
bath model, and the ring coupled to a spin bath, is that 
in the oscillator bath system, decoherence is in a certain 
sense a mere side effect of the dissipation taking place 
each time the particle excites an oscillator. On the other 
hand in the spin bath model, no such dissipation occurs, 
only phase decoherence. This difference is most obvi- 
ously seen in the centre of mass dynamics of a particle 
wave-packet - for the oscillator bath model a wave-packet 
initially moving around the ring will dissipate centre of 
mass momentum (formally this happens via the interac- 
tions between qr and cr^' in the influence function) , slowly 
bringing it to rest. However, as we see in the next sec- 
tion, for a ring particle coupled to a spin bath, the centre 
of mass momentum of a wave-packet is completely con- 
served, even in the strong decoherence limit, provided 
the spin bath dynamics is governed by its coupling to 
the ring particle (the typical case). This leads to some 
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counter-intuitive features, as we now see. 



V. WAVE-PACKET INTERFERENCE 

It is interesting to now turn to the situation where two 
signals are launched at t = from 2 different points in 
the ring. The idea is to see how the spin bath affects 
their mutual interference, and how, by effectively cou- 
pling to the momentum of the particle, it destroys the 
coherence between states with different momenta. We 
do not give complete results here, but only enough to 
show how things work. 

We therefore start with two-wave-packets which will 
initially be in a pure state, and will then gradually be 
dephased by the bath. In the absence of a bath, we will 
assume the wave function of this state to be the symmet- 
ric superposition 



(68) 



where we assume the usual symmetric ring with flux 

and Z — \Jj2n=o 6^''"'"^'^/^)^-° is the wave-function nor- 
malization factor. At t — 0, one of the packets is centred 
at the origin, and the other at site jo, and they both have 
width D. Note that the velocity of each wave-packet is 
conserved, and at times such that Aot = 2n, they cross 
each other. From (IMj) we see that the main effect of 
the flux is to shift the relative momentum of the wave- 
packets. It also affects the rate at which the wave-packets 
disperse in real space - this dispersion rate is at a mini- 
mum when ^/N — ^. 

The free-particle wave function in real space is then 



1 



N-l 



X (e' 



-(fe„-7r/2)^£)/2 



-2iAat cos {k„+i/N) 



(71) 



g-'ijkn g-2iAot cos {k„-i/N) 



where the two wave-packets are assumed to have Gaus- 
sian form: 

N-l 

\Mt)) = I E e-C^"--/^)^^/^ 

n=0 

X g-»Jofc„-*2Aotcos(fc„-^<I./JV)|^^^ (gg) 



SO that the probability to find a particle at time t on site 
J is P(j) = \^,{t)\\ ' 

Let us now consider the effect of phase decoherence 
from the spin bath. Using the results for Pjj' {t) from the 
last section, with an initial reduced density matrix 



-(fe„-ir/2)^D/2 



n=0 



X e-*''^°*™"('="-*/^)|27r-fc„) (70) 



(in) 



l*j(i = 0))(*,'(t = 0)| 



(72) 



we find a rather lengthy result for the probability that 
the site j is occupied at time t: 



TV— 1 +CX; 

^^■W-^ E E e-«'='.-/2)^+('="'"-/^)^)^/^^^o(™) 

rt,n'=0 m= — CO 

x{e*0--Jo)(fe„-fc„0 j„^(4A„fsin((/fc„ - fc„,)/2))e"""''"+''"'^/'+*/'^^+ (^3) 
_^g-»(fe„-fe„,b j^^(4A„isin((fc„ - fc„,)/2))e™(^'="+'="'^/'"*/^^ + 
+ [e'«J'-^")'="+^'="'' J„(4Ao<sin((fc„ + A:„0/2))e""((^"-'^"'^-*/^' -f H.c.]} 

One can also, in the same way, derive results for the current in the situation where we start with 2 wave-packets. 
We see that expressions like ((73)) are too unwieldy for simple analysis. However in the strong decoherence limit ((73)) 
simplifies to: 

P^it) = ^ e-«'=--/2)^+('="'--/2)')^/2{e^W-^o)(^-"-'='^')Jo(4Aisin((fc„-fc„,)/2))+ 

+ e-'J^'^"-*^"') Jo(4Atsin((fc„ - + [e'^^i-3o)k^+3K-) j^{^4^/s,t sin ((A:„ + K,)/2)) + H.c]} (74) 

I 



and again we see that the flux has disappeared from this equation. 
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FIG. 6: Colour online: Plot for Pj{t) as a function of both j 
and Aot in the strong decoherence limit, jo ~ 50 and — 
100. The relative velocity is -I, in phase units. Top: global 
view. Bottom: a particular peak 



This result is shown in Fig. |6l As one might ex- 
pect, the interference between the two wave-packets is 
completely washed out in this strong decoherence limit. 
However there is also a more unexpected feature - each 
wave-packet now has portions moving in opposite direc- 
tions to each other. The explanation is to be found by 
noticing that as the particle hops, at the same time caus- 
ing the bath spins to make transitions, the topological 
phase it exchanges with the spins also changes the to- 
tal phase around the ring seen by the particle. Thus, 
from the point of view of the particle, these transitions 
are forcing the total flux through the ring to fluctuate, 
in a way which depends on the trajectory followed by 
the particle. This dependence, such that the changing 
phase is conditional on the particle path, is of course 
why we get decoherence. Now the changing effective flux 
changes the particle momentum and velocity, and in the 
case of the pair of wave-packets here, it also changes their 
relative momentum. Indeed, given that the transforma- 
tion $ — > $ + TT completely reverses the momentum, we 
see that a strong coupling to the bath spins can even 
cause a part of the initial wave-packets to reverse its di- 
rection. However we emphasize that the centre of mass 
momentum for the combined wave-packet system has not 
changed - the average momentum imparted to the ring 
particle is zero. Thus, as noted in the last section, the 
decoherence caused by the bath is not accompanied by 
any net dissipation of the particle momentum, or of its 
energy. Indeed, if a wave-packet starts off with a net an- 
gular momentum around the ring, this will be conserved. 



long after all coherence has been lost. 

Note that these results are not the same as one would 
get by just adding a fluctuating noise S^{t) to the static 
flux. Such an external noise term will also cause 'noise' 
decoherence, but of a quite different form from that of the 
spin bath decoherence discussed here, since there is no 
correlation between the noise and the particle dynamics. 
In fact, as we will discuss elsewhere, a fluctuating flux 
noise acting on the particle causes exponential decay in 
time of the particle correlation functions, quite different 
from the power law decay typical for the present case. 



VI. SUMMARY & CONCLUSIONS 

Let us flrst recall the main results derived in sections 
II-V. In section II we show how the basic Hamiltonian ([1]) 
we have studied can be derived, including the fairly severe 
approximations that are involved. No attempt is made 
to connect the model with any speciflc physical system, 
since the main focus of this paper is to study pure phase 
decoherence in a solvable model. The most important 
features of the model are that (i) the spin bath which 
couples to the model can cause severe phase decoherence 
with no dissipation, and (ii) the phase interference in the 
ring (including Aharonov-Bohm oscillations) is affected 
in a rather fascinating way by the decoherence; and (iii) 
the model can be solved exactly. The main tasks we set 
ourselves in this paper were to set up a formal apparatus 
to solve this model, and to study some aspects of the 
decoherence dynamics in it. 

Before studying the decoherence, it turns out to be 
important to develop the detailed solution for the dy- 
namics of the iV-site ring without the bath, in section 
III - to our surprise, this does not seem to have been 
done before. It is convenient to develop the free particle 
density matrix p°{t) and its propagator K°{t) as double 
sums over winding numbers around the ring (see eqtn. 
(|A9[) for the propagator); but we also show how this can 
be rewritten as a single sum over winding numbers (see 
eqtn. (jAlOl) . a form more useful for numerical work on 
large rings. 

The importance of the work on the free particle prob- 
lem is seen in the result (|^ for the propagator IC{t) 
of the reduced density matrix once one integrates out 
the spin bath - we can find this by summing over wind- 
ing numbers an expression involving matrix elements of 
K°{t) and a weighting function Fjj/(p,_p'), the influence 
function. The influence function can be found exactly 
(Appendix A. 2); we do this both for the case where the 
couplings between the particle and the bath spins are 
flxed, and the case where we make an ensemble average 
over these couplings. This allows us to derive a whole 
series of exact expressions for the propagator ]C{t) (see 
Appendix A. 2, eqtns. (|A26|) and (|A30p ). and thence for 
the time evolution of the reduced density matrix, the 
probability density, and the current (section IV). In de- 
riving results for these physical quantities, one finds that 
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the full details of Fjj,{p,p') are not required, but only 
certain matrix elements (often only the element Fq(p), 
produced by letting /i = j — j' + 1 — I' = 0, and p' —p = p; 
see (|44|) et seq.). The characteristics of the decoherence 
are controlled by these. 

It turns out that the decoherence dynamics, and how it 
affects different physical quantities, depends on the ring 
size N, the flux $, and the form of the function Fq{p); and 
also on the initial state of the system. A full exploration 
of this large parameter space would take a lot of space, 
so we have focussed on certain questions. One of these 
is the flux dependence of the various physical quantities, 
and how phase decoherence affects these - how this works 
is shown in section IV, with details given for a 3-site ring. 
We also look at the interference of 2 wave-packets on a 
large ring, in section V. This of course depends crucially 
on the flux, and as we switch on decoherence, this flux de- 
pendence disappears, even though the wave-packets still 
propagate (although the coupling to the spin bath also 
strongly distorts the shape of the wave-packets). In all 
cases we find that the detailed dynamics is quite different 
from what one would get if the decoherence was simu- 
lated by adding flux noise to the problem - in particular, 
all coherence properties show power law decay in time, 
instead of exponential decay. 

Because of the size of the parameter space, there is 
much that is unexplored in this paper - in particular, we 
expect the decoherence dynamics, and its dependence on 
flux, to depend very dramatically on the form of Fq{p)] 
and we have hardly explored the dependence on ring size. 
Nor have we attempted any connection to experiment. 
The main reason for this is that for a detailed compar- 
ison with experiments on most real systems one has to 
add two crucial ingredients, viz. (i) we must add back 
the local flelds {hfe} acting on the bath spins, and the di- 
agonal couplings {7;^} (compare eqtns. and ([7])); and 
(ii) in many physical applications there are also impor- 
tant couplings to delocalised modes like phonons, which 
are modeled using oscillator bath interactions. Actually 
one can also solve the problem when these couplings are 
added, in certain parameter ranges - this will be the sub- 
ject of future papersi^ 
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text, and also explain some of the mathematical trans- 
formations required to go from single sums over winding 
number to double sums. 



1. Free Particle 

We consider flrst the free particle for the A^-site sym- 
metric ring, with Hamiltonian 



(Al) 



<tj> 



and band dispersion = 2AoCos(fc„ — ^/N). 

For this free particle the dynamics is entirely described 
in terms of the bare 1-particle Green function 



q,,(i)^(j|G°(i)|/)^(j| 



N 



-i2Agtcos(k^~<S>/N) ik^{j'-j) 



(A2) 



which gives the amplitude for the particle to propagate 
from site j' at time zero to site j at time t. This can be 
written as a sum over winding numbers m, viz.. 



( i^ptY i'S>/N(l-2m) 



e=o rn=a ^ ' 



-E 



(A3) 



71=0 



This sum may be evaluated in various forms, the most 
useful being in terms of Bessel functions: 



N-l +00 



^E E Jm{2Aot){-z) 



771 im{kn — ^/N)-\-ikn{j—j ) 



n—0 m— — 00 
+00 



J2 Jm i2Aot) (-i )™e-™*/^(5^p,^+,_ 



m— — 00 
+c>o 



-i(7Vp+i'-i)(*/iV+7r/2) 



P— — 00 



(A4) 



(this last form, where we have eliminated the sum over 
winding numbers, is also of course directly derivable from 
^2[) ). We can also write this last form as 

Gn' it) = E e^^*+*^(^-^"')/7Vp+,-,' h2^AJ) , (A5) 



Appendix A 

In this Appendix we derive some of the expressions for 
Green functions and density matrices that are used in the 



where we use the hyperbolic Bessel function /q(x), de- 
fined as Ia{x) = {i)~°'J^{ix). 

Consider now the free particle density matrix. As dis- 
cussed in the main text, we have in general some initial 
density matrix p^™' = {l\p°{t = 0)\l') at time t = 
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(where I and I' are site indices). Then at a later time t 
we have 

Pjj'it)^J2^n'M'WfP (A6) 

where K°-, (t) is the propagator for the free particle 
density matrix. Its form follows directly from the defini- 
tion of this density matrix as p°{t) = \^(J{t)){■^|J{t)\, where 
\ip{t)) is the particle state vector at time t. One thus has 

^j;-',n'W = q/Wq'/'W^- (A7) 



An obvious way of writing this propagator is then: 



pp' 



where we have a double sum over winding numbers p,p'. 
The explicit form for K°^, n,(p,p'\t) is then given from 
(Ei) as 



pp' 

= ^ e»(f-f')*e^*0-^-'+'-O/A^/^^^^._^(_2zA„t)/^p,+^.,_,(2zA„t) . (A9) 
pp' 

However this expression is somewhat unwieldy, particularly for numerical evaluation, because of the sum over pairs of 
Bessel functions. It is then useful to notice that we can also derive the answer as a single sum over winding numbers, 
as follows: 

1 

K° \^ p-«(fe„O'-0-fe„/O ))+4iA„tsin[$/Ar-(fc„+fe )/2]sin[(fe„-fe^,)/2] 

n.7i'—0 
^ N-1 00 

= E E ^p[4A„isin(fc„/2)]e^f(*/^-'='^+'="./2)-^'=.(^-')+*('=.-'=-.)(^'-'') 

n,m— P— — 00 
00 N-1 

E Y.i'W+r-:+i-i'i^^otsHkJ2)y'^^^(^+^'-^-^^^ . (Aio) 

p' — — oo m—0 

I 



In the second step we replaced n' ^ m — n. In the 
third step we also used the identity J2n'=o e*'^-'^ = 

The result for the density matrix then depends on what 
is the initial density matrix, according to (jA6|) . If we 
start with p*^™) = |0)(0|, the density matrix is then just 
Pjj'^^^ ~ ^jj'.ooi^)- then gets a much simpler ex- 

pression; the density matrix at time t is: 

^ N-1 oc 

'^jVW = ]^E E •^iVp'+/-,[4A„isin(fc„/2)] 

rn—O p' — — oo 

X f.iHNp' +j' -3)-ik^U+j' -Np')/2 

(All) 

It is useful and important to show that the double- and 
single-sum expressions (|A9P and (jA10|) are equivalent to 
each other. To do this we use Graf's summation theorem 
for Bessel functions, in the form: 

+00 
J,(2a;sin-)(-e-*«)5 = ^ J,+^(x) J^(x)e*^« (A12) 

/J,— — 00 



We set 9 = 0, , ...gf^, ... which is the in 

(jAlOp and multiply by e"**-' on each side. We then have 

j^(2xsin^)e-^('="+^'te-"="^' 

= ^.+,.(a^)^A<(a;)e^^^"^"™ (A13) 



Noticing then that X]m=o e**^"" = N J2p ^Np.n we then 
do the sum over m; only — j — Np survives, and thus 

^ ^ J,(2a;sin^)e~'('='"+'^)5e-^'='"^' 

m=0 

= ;^ E Jnp+j+^{x)Jnp+j{x) (A14) 
p 

Setting V = Np' + j' - Np- j + 1 - I' , x = 2 A^i, we then 
substitute back into (|A9p . to get 



16 



pp' 

= ^ E e+^(^^+^'-^+'-'')( t+f ) J^,+,.,_,-+,_, (4A„t sin 

p m—O 
p m— 

The propagator p is Hermitian, ie., K°j, = iir",^- iri{t)*; setting p' = —p, we then have 

K^M'i^ - ^E E J-^p'+,-/+.-K4A„tsin^)e^(^^'+^'-+'-'')f +..™0-+/-.-.'-Arp')/2 

p' m—O 
7V-1 00 

= ^ E E ^^p'+.'-,+.-.[4Ao^Bin(fc,„/2)]e^^(~^'+^'-^+'-'')-^„0-+/-^-r-^p')/2 (^16) 

m—O p'^ — cxD 

where in the last hne, we set k„i — > —km, and use the fact that for integer order n, J„(— x) = J_„(a;). Thus we have 
demonstrated the equivalence of the single and double sum forms for the density matrix. 



2. Including Phase Decoherence 

To calculate the reduced density matrix for the par- 
ticle in the presence of the spin bath, we need to av- 
erage over the spin bath degrees of freedom. We will 
do this in a path integral technique, adapting the usual 
Feynman-Vernonl^ theory for oscillator baths to a spin 
bath; the following is a generalization of the method dis- 
cussed previously.lS We can parametrize a path for the 
angular coordinate Q{t) which includes m transitions be- 
tween sites in the form 

m 

Q(m){{q^},t) = Q{t = Q) + Y,q,e{t-U) (A17) 

i=l 

where 0{x) is the step- function; we have transitions cither 
clockwise (with qj = +\) ov anticlockwise (with qj = — 1) 
at times ti,t2, ■ ■ ■ ,tm- The propagator /C(l,2) for the 
particle reduced density matrix between times ri and T2 
is then 

/c(i,2) = r^de r'de'e-*(^°[«i-«°[«'i)j-[e,e'] 

(A18) 

where S'o[6] is the free particle action, and ^"[9, 9'] is 
the "influence functional" defined by 

J-[9,9'] - ll{Uk{e,t)ul{Q\t)) (A19) 

k 

Here the unitary operator Uk{&,t) describes the evo- 
lution of the k-th environmental mode, given that the 
central system follows the path Q{t) on its "outward" 
voyage, and Q'{t) on its "return" voyage. Thus J"[9,9'] 
acts as a weighting function, over different possible paths 



(0(t), Q'{t)). The average (...) is performed over environ- 
mental modes - its form depends on what constraints we 
apply to the initial full density matrix. In what follows 
we will assume an initial product state for the full parti- 
cle/environment density matrix. 

For the general Hamiltonian in eqtns. ©-(U), the en- 
vironmental average is a generalisation of the form that 
appears—'—, when we average over a spin bath for a cen- 
tral 2-level system, or "qubit". The essential result is 
that we can calculate the reduced density matrix for a 
central system by performing a set of averages over the 
bare density matrix. For a spin bath these can be re- 
duced to phase averages and energy averages; and for 
the present case it reduces to a simple phase average. 

To see this, notice that the sole effect of the pure phase 
coupling to the spin bath is simply to accumulate an ad- 
ditional phase in the path integral each time the particle 
hops. Just as for the free particle, we can then classify 
the paths by their winding number; for a path with wind- 
ing number p which starts at site I (the initial state) and 
ends at site j, the additional phase factor can then be 
written as 

/ (NO) (j-lj> \ 

exp{-zp^ E - E ("r-'T.)} 

*: \(mn) = {01> {mn) = {l,l+l) J 

(A20) 

and for fixed initial and final sites, this additional phase 
only depends on the winding number. 

Consider now the form this implies for the reduced 
density matrix of the particle, once the bath has been 
averaged out. The equation of motion for the reduced 
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density matrix will be written as 
1,1' 

where JCjji^ii'{t) is the propagator for the reduced density 
matrix. Now the key result is that 

fCn'M'it) - E (P^P'-^ ^)^"'' (P^P''^ 

pp' 

where the function K°j, ii,{p,p';t) is the free particle 
propagator for fixed winding numbers p,p' (see eqtn. 
(jASp above). This form follows from the argument just 

I 



given, viz., that the only effect of the spin bath is to add 
the extra phase factor (jA20|) in each path in the path in- 
tegral for the propagator. Thus the influence functional, 
initially over the entire pair of paths for the reduced den- 
sity matrix, has now reduced to the much simpler weight- 
ing function Fj', {p,p'), which we will henceforth call the 
" influence function" . To evaluate this influence function, 
we must specify what kind of bath average we wish to 
take. We consider here the two cases discussed in the 
text, viz., (i) a simple average < ... > over bath states, 
and (ii) an average << .... >> over both the bath states 
and over a distribution P(a™") of couplings to the bath 
spins. The results are obtained as follows: 



(i) In the case of fixed bath couplings a™", the average is obtained by simply inserting the phase factors ()A20p 
from above into the paths of different winding number. We then get: 

In the symmetric coupling case where ct™" a^, this expression reduces to a much simpler result: 

Ffy{p,p') = (e*(M+^p)i:.a.-T.^ (A24) 

where we define 

fi = f-j + l'-l, p = p'-p. (A25) 

If the particle is launched from the origin, this gives the even simpler result (|48p quoted in the main text. 
We can now give the explicit result for the propagator of the reduced density matrix in double sum form as 

= Y.ie^^^+^P^^^^"-^") e-P)*e-'*'^/^z^P+^J^^+^._,(2A„i)J^^,^^.,_,,(2A„i) (A26) 
pp' 

Typically the average < ... > here will be over a set of thermally weighted states, but this is not required - in principle 
one could average with a non-thermal state of the bath (or even a definite bath state, eg., one which had been polarized 
beforehand - in this case no bath averaging is required at all). If we do assume a thermal state, then since all bath 
states are degenerate, and hence equally populated at any finite T, eqtn. (jA24p reduces to 

F"-'{p,p') - {e^(>^+Np)i:.a.-o:) ^ l[cos{{Np + fi)\cx\) (A27) 

k 

These results can also be written in single sum form - we do not go through the details here. Both forms are fairly 
easily summed numerically, even for rather large rings. 

(ii) In the case where we must also average over the distribution of bath couplings, the same phase factors appear, 
but now we must average over the bath spin couplings; this gives instead 

F^Up^p') = n fdc^rPic^D 

k 

X (e-^(p-p')ES:„-,^iTo> «r"-T.g-«E<^'2<.,.+,> «r -T. "r-^.^ ^^28) 

where we put a bar over the infiuence function to signify an extra average over bath couplings; again < ... > signifies 
the average over bath states, and P(q;™") is the probability weighting for the different bath couplings. Typically 
it makes little sense, in this ensemble average, to have any dependence of the P(q;™") on the link mn, so that this 
reduces to 

Ff;,{p,p') = ]^((e«(A'+A'p)a.-'T.)) = [] / dakPiau) (g^Ca'+^'p)"-'^^') (A29) 
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where H = j' — j + I' — I, and p = p' — p as above, and where we use the fact that the distribution P{ak), defined 
for the individual bath couphngs, reduces in an ensemble average to a simple average over some coupling strength a, 
with weighting P{(y.), acting on some representative spin tr (see main text, section IV. B). Thus our explicit result for 
the propagator of the reduced density matrix is now 

pp' 

To make this result more specific, let us assume the Gaussian distribution of couplings given in the text (eqtn. 
()45p ). If we again assume a thermal average then we can easily evaluate this average, in the same way as in the main 
text, to get 

{{e'^t^+ma-(T^^N^ ^ cxp[-A(iVp + Ai)V2] (A31) 

As before, the result (jA30p can be rewritten as a single sum, and this is again fairly easily summed numerically. We 
note that the effect of the influence function is now to rapidly suppress paths in which \Np + fi\ is non-zero. 

I 



Now consider the current Ijj^i{t). Again we may dis- 
tinguish between the case where we average only over 
the bath states, and the case where we average over both 
the bath states and the bath couplings. In the case of a 
single average over bath states, with fixed couplings, the 
current is given by Eqs. ((57)) and (l58t . with the brack- 
ets in (|57p replaced by double brackets when we average 
over the bath couplings as well. We see that formally ev- 
erything depends only on the phase between sites j and 

I 



j + 1, via the bath-generated phase-dependent coupling 
Ajj+i, and on the density matrix element pjj^i{t) and 
its conjugate at time t. However this apparent simplicity 
is deceptive, because the density matrix depends itself on 
the form of the initial density matrix at time t = 0, 
and on the propagation of this density matrix in the in- 
terim; thus ([57]) contains implicitly the full propagator 

kIUp,p'). 



Using the results derived above for this propagator, we can now derive expressions for Ijj+i{t). In what follows we 
only quote the results of the case of a bath average with fixed couplings - the case where one also does an ensemble 
average over the couplings is easily deduced from these expressions, following the same manoeuvres as above. The 
results can be found in both single and double winding number forms. The double Bessel function form is 

2Ao ^ J7Vp+j-i(2Aoi) Jjvp'+j+i_('(2Aoi) 
pp' 

X Re(p|;^)^^(f-f')e'[(P-P')+^l* e^^^'^'^^" ^Z^-m^ "^"■'^.^2.1:. E<r„>4„,'+i> (A32) 

Again, let us make the assumption of a completely ring-symmetric bath, so that cy.^^ — otk- Then we get 

= 2A„ ^ ^ JNp+,-ii2Aj)JNp'+j+i-i'i2Aj)Fi^i,ip',p) x Re[p[;r^e^*[p'-f+('-'')/^)l] (ASS) 
pp' 1,1' 

From this we can derive the single Bessel Function summation form as follows. Using the equation 

Y.JNp+n-i{x)JNp+n-i+u{x) = ^ E Jfe (2:E siu ^)e-^("-')'='" -^C^" "^^''/^ (AS4) 

p ni—0 

which is another form of Graf's identity we set v = N{p' — p) + l + l~l',x = 2Aot; then 
2Ao 

Ijj+iit) = -j;f-'^'^'^JNp+i+i~i'{'i^otsiii-^) 

m=0 p 1,1' 

X e~^'='"[^+"-('+'')/2l^JVp+i+'-''Fj,,(p)Re[p|;r^e'*[(p'-f+'-'')/^)l] (AS5) 

where we define Fui {p, 0) = Fw (p). 
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If we make the assumption that the particle starts at the origin, these results simplify considerably; one gets 

i2Aj)Foip',p) cos[{-N + <^){p' - p)] 

pp' 

- E E•^^P+l(4Ao^sin^)e-^-(^+^■)^^^+lFo(p)cos[(|iV + $>] (A36) 

m— p 

I 

for the double and single sums over winding numbers, is used in the text for practical analysis, 
respectively; and Fo {p) = Fjj (p, 0) . The latter expression 
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